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RESEARCH MEMORANDUM 


LATERAL AND DIRECTIONAL DYNAMIC -RESPONSE CHARACTERISTICS 
OF A 35° SWEPT -WING AIRPLANE AS DETERMINED 
FROM FLIGHT MEASUREMENTS 
By William C. Triplett and Stuart C. Brown 


SUMMARY 


Lateral and directional dynamic -response characteristics of a 35° 
swept-wing fighter-type airplane determined from flight measurements are 
presented and compared with predictions based on theoretical studies and 
wind-tunnel data. Flights were made at altitudes of 10,000 and 35*000 
feet covering the Mach number ranges of 0.50 to 0.8l and 0.50 to 1.04, 
respectively. Recorded data consisted of transient responses in yawing 
velocity, rolling velocity, and sideslip angle to pulse-type motions of 
the rudder and of the ailerons. These transient data were converted 
into frequency-response form by means of the Fourier transformation and 
compared with predicted responses calculated from the basic equations of 
motion. The equations, or transfer functions, that best describe the 
various measured responses were evaluated by a curve-fitting process 
involving the use of templates and an analogue computer. By this method 
it was generally possible to find equations, of simple form, that closely 
matched the experimental frequency responses between 1 and 10 radians 
per second and at the same time adequately described the recorded time 
histories . 

Experimentally determined transfer functions were used for the 
evaluation of the stability derivatives that have the greatest effect 
on the dynamic response of the airplane. The values of these derivatives, 
in general, agreed favorably with predictions over the Mach number range 
of the test. There were notable exceptions, however, in some cases at 
the flight altitude of 10,000 feet. These discrepancies are attributed 
to aeroelastic deformations of the wing and tail. 


Another departure from theory was disclosed 
spiral root which was many times larger than expecft 
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INTRODUCTION 


This report describes the second phase of a flight-test program 
being conducted by the NACA for the purpose of determining the dynamic - 
response characteristics of a 35° swept-wing fighter-type airplane. The 
first phase of this program, which considered only the longitudinal case 
was reported in reference 1. The present report applies the methods of 
reference 1 to the lateral- and directional -response characteristics of 
the test airplane. Frequency responses and the associated aerodynamic 
derivatives are evaluated from records taken throughout the Mach number 
ranges of 0*50 to 1.0^ at an altitude of 35*000 feet and 0.50 to 0.8l 
at 10,000 feet. 

The test procedures and analysis methods used are essentially the 
same as described in reference 1. Transient responses to both rudder 
and aileron disturbances are measured and analyzed to give frequency 
responses of yawing velocity, rolling velocity, and sideslip angle. 
Responses to transient rather than sinusoidal control inputs have been 
chosen for analysis because of convenience in making flight measurements 
While certain aerodynamic information can be determined directly from 
transient time histories, the use of the frequency response concept 
allows a more complete evaluation of the dynamic behavior of the air- 
craft. The effects of different modes on the over-all motion can often 
be shown more clearly in the frequency plane, especially when these 
motions are complicated by structural deformation. 

In this investigation, transfer functions with numerical coeffi- 
cients are obtained directly from the measured frequency responses. 
Supplementary calculations are then made to evaluate those aerodynamic 
derivatives that exert the strongest influence on the airplane response. 
Data of this type are of particular interest in the study of airplane - 
autopilot .combinations. Knowledge of the airplane transfer functions is 
necessary in the determination of the dynamic characteristics that are 
required of an autopilot to satisfactorily control the aircraft. 

Wherever possible the results of these tests are compared with 
predictions based on wind-tunnel and theoretical data and also with the 
results of other flight tests. 


NOTATION 


0^ lift coefficient 

Cl rolling- moment coefficient 
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C n yawing -moment coefficient 
Cy side— force coefficient 

D the operator 

I x moment of inertia about the X axis, slug-feet squared 

Iy moment of inertia about the Z axis, slug-feet squared 

IxZ P ro <i u ct of inertia, slug-feet squared 

R,I real and imaginary parts of a complex quantity 

5 wing area, square feet 

V velocity, feet per second 

W weight of airplane, pounds 

b wing span, feet 

g acceleration due to gravity, feet per second squared 

i J-T 

m mass of airplane, slugs 

p rolling velocity, radians per second 

q dynamic pressure, pounds per square foot 

r yawing velocity, radians per second 

t time, seconds 

a angle of attack, degrees 

3 sideslip angle, radians (except as noted) 

7 flight path angle, degrees 

6 control deflection, radians (except as noted) 

6 a total aileron deflection, radians 
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k 

6 r rudder deflection, radians (except as noted) 

£ damping ratio 

A, root of the characteristic equation 
a real part of a complex root 

9 angle of hank, radians 

0 phase angle, degrees 

•\|r angle of yaw, radians 

03 frequency, radians per second 

U) 1 natural frequency of oscillation, radians per second 
0 ) n undamped natural frequency, radians per second 

% 5W2VT’ Ier radlan 


SC Z 


C lr 


per radian 

% 

dc 7 

ST' per 

radian 

C °P 

dC n 

per radian 

d(pb/2V)» 

Cn r 

ac n 

per radian 

d(rb/2v)» 

c n p 

acn 

ST’ per 

radian 

% 

ST' per 

radian 

% 

dCy 

per 

radian 

c n 6 

ac n 

per 

radian 
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C *5 

, per radian 

c 2 . 

ip 

ac 2 

, seconds squared per radian 

Lp 

2Vfy C 2 p’ per Sec ° nd 

L r 

2VI X C *r> per second 

• L P 

tj— ^ P er secon( i squared 

% 

qSb 2 

Cnp, per second 

N r 

gSb 2 , 

Cn r > P er second 

Ng 

Cnp> P er second squared 

Y P 

Sf C Yp> P er second 

p 6 

C 7 y per second squared 
-*-X S 

N 8 

Cng> P er second squared 

Y 5 

^ Cv , per second 
mV * 5 * 

r X 

I XZ 

ix 

r z 

x xz 

*z 

K 1 

£ cos 7y per second 

• 

k 2 

^ sin 7 , per second 
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TEST EQUIPMENT 


The test airplane was a standard North American F-86A-5 with exter- 
nal instrument booms added as shown in figures 1 and 2. The physical 
characteristics of this airplane are described in table I. 

Standard NACA instruments were used to record airspeed, altitude, 
rolling and yawing velocities, normal acceleration, angle of attack, 
sideslip angle, and rudder and aileron positions. All recordings were 
synchronized at 0.1 -second intervals by a common timing circuit. The 
true Mach number was obtained from the nose -boom airspeed system 
described in reference 2. 


Rate gyros were used to measure yawing and rolling velocities about 
the reference axes of the airplane. The yaw rate gyro had a range of 
±0.5 radians per second and a natural frequency of 10 cycles per second. 
Corresponding values for the roll rate gyro were ±2.0 radians per second 
and 18 cycles per second, respectively. In both cases the damping ratios 
were approximately O.J. Sideslip angles were measured by a vane-type 
pickup and recorded on an oscillograph. Rudder and aileron deflections 
were measured by control position recorders that were linked directly 
to the control surf aces . 

The dynamic characteristics of recording instruments are of extreme 
importance in investigations of this type. A flight record in general 
contains the combined response of the airplane and recording instrument. 
If the instrument has a linear second-order response with known damping 
ratio and natural frequency, then its response can be subtracted from 
the combined response in the frequency plane. Obviously, it is desira- 
ble to use instruments with characteristics such that the necessary 
corrections are a minimum (i.e., high natural frequency and damping ratio 
of approximately 0.7). 
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The two rate gyros (rolling and yawing velocities) used in the 
present investigation were considered to he very satisfactory in this 
respect and no corrections were applied to the data. Because the 
response of the sideslip vane, however, was unknown and suspected to he 
nonlinear, results obtained from this instrument were not considered 
reliable at high frequencies. The control position recorders (aileron 
and rudder) also had unknown frequency responses hut tests of similar 
installations have indicated very high natural frequencies, so they 
were assumed to give valid records over the frequency range of interest. 


TEST PROCEDURE 


The flight procedures consisted essentially of recording airplane 
responses to both aileron and rudder disturbances. Pulse-type inputs 
were used in all but one flight run. All flights were made at altitudes 
of 10,000 and 35>000 feet in the Mach number ranges of 0.50 to 0.8l and 
0.50 to 1.0^-, respectively. The corresponding trim lift coefficients 
varied from 0.17 to 0.07 and 0.51 to 0.12, respectively. 

When one control was specified as the disturbing element, the other 
control was held fixed during the entire maneuver. After application of 
the pulse input, both controls were held fixed until the oscillatory 
motion of the aircraft had essentially subsided. The airplane responses 
used in the analysis were rolling velocity, yawing velocity, and side- 
slip angle. Sample time histories of the responses of the quantities 
to the appropriate control inputs can be seen in figure 3* 

All flight runs at speeds below a Mach number of 0.95 were made in 
trimmed level flight, but to obtain data at the higher speeds it was 
necessary to dive the airplane. Flight altitudes changed as much as 
2,000 feet during each diving run although there was little variation 
in Mach number. For analysis purposes the altitude and dynamic pressure 
were assumed to be constant at their average values during each run. 


METHOD OF ANALYZING EXPERIMENTAL DATA 


The basic method used in the analysis of the transient flight data 
was described fully in reference 1 but is briefly summarized here. 
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The evaluation of the transfer functions that describe the motions 
of the airplane were carried out as indicated in the sketch: 



The transient responses were first corrected wherever necessary so 
that the measurements would conform to the stability axes notation that 
is commonly used. This correction is explained in appendix A. The data 
were then analyzed to obtain the Fourier transformations by means of an 
IBM machine calculating method employing an adaption of Simpson 1 s rule. 
These calculations were carried out at a number of frequencies between 
0. 5 and 1 6 radians per second. The resultant frequency responses were 
plotted as shown in the form of amplitude ratios and phase angles. 

As discussed in appendix B, the type of input used in flight defi- 
nitely places a limit on the accuracy of the Fourier analysis of a tran- 
sient record. In general, to obtain the widest usable frequency range 
a pulse-type input should be used. When low frequencies (below 1 radian 
per second) are desired, a step input is preferable although this type 
of disturbance may result in motions that exceed the ranges of linearity. 
The method used, while subject to these limitations, is extremely accurate 
provided that the flight records are tabulated at enough time intervals 
to clearly define the data. 

The second step consisted of fitting the graphical frequency 
responses to a set of dynamic response templates which define the first- 
and second-order complex functions 1 + iu and 1 + 2^iu - u 2 where u is 
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a nondimens ional frequency variable and £ the damping ratio. This set 
of templates , as described in reference 3* contained one pair of first- 
order curves (amplitude and phase) and a family of second-order curves 
with different values of the parameter £ . By combining templates in 
the proper manner it was generally possible to find one combination that 
would closely match both amplitude ratio and phase-angle plots. This 
particular combination defines the equation, with approximate numerical 
coefficients, that most nearly describes the graphical frequency response. 
This equation is normally referred to as the transfer function. 

The final step of the analysis involves the use of an analogue com- 
puter on which the airplane transfer functions are set up. If a time 
history of an actual recorded control motion is supplied as an input to 
the computer, then the output should be indentical to the response 
measured in flight. Since the templates give only approximate values 
of the coefficients, these values were altered until the output of the 
computer satisfactorily matched the experimental time response. Thus, 
in addition to refining the results, the computer automatically furnishes 
a check on the previous computations . 

After the transfer functions that describe the various airplane 
responses were obtained, the stability derivatives C nR , C Zr , CV, C z , 
c Z &a > and C n were evaluated as explained in a later 1 section. ^ P 


RESULTS AND DISCUSSION 


The results discussed in the following paragraphs were obtained 
from transient time histories as explained in the previous section. 
Flight evaluated frequency responses, transfer functions, and stability 
derivatives are presented and compared with predictions based on wind- 
tunnel data and theoretical studies. 


Frequency Responses 


Plotted in figure V are typical flight evaluated frequency responses 
of rolling velocity, yawing velocity, and sideslip angle to rudder anri 
also to aileron inputs. These were all obtained at an altitude of 
35*000 feet; responses for 10,000 feet showed similar characteristics 
and have not been plotted. The purpose of these figures is to show 
general trends with varying Mach number , and therefore smooth curves 
have been faired through the calculated test points. In most cases more 
than one flight record was analyzed at each flight speed in order to 
check the data for consistency. Only at the highest test speeds (above 
a Mach number of 1.0) was there appreciable inconsistency. 
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The responses to rudder inputs have not been plotted at frequencies 
greater than 8 radians per second because there was considerable scatter 
and also a lack of well-defined trends in the data at the higher fre- 
quencies. The aileron responses, however, are shown to 16 radians per 
second. Wherever necessary for clarity or because of erratic data, 
parts of some of the curves have been omitted. The p/b a response is 
shown at only three speeds because of a failure in the sideslip-angle 
recording system. 

With minor exceptions, the curves show consistent and gradual vari- 
ations with Mach number. One such exception can be seen in figures 4(d) 
and 4(e) where the amplitudes of p/b a and r/S a at a Mach number of 0.6l 
lack the customary resonant peaks. This is the result of time histories 
in which there was no oscillatory motion. This unusual characteristic 
can be explained by reference to the predicted transfer functions devel- 
oped in appendix C and discussed in the following section. The predicted 
p/5 a response equation for a Mach number of 0.6 at 35,000 feet is 

p_ = 21 . 1 ( D^+0 . 455D+6 . 74 ) 

S a ( D+0 . 00113 ) ( D+2 . 20 3 ) ( D 2 +0 . 438 D+ 7 • 25 ) 

It can be seen that the two quadratic terms are nearly identical and 
thus the oscillatory mode is effectively canceled. Since rolling and 
yawing motions are coupled, the r/6 a response must exhibit the same 
characteristics at this particular speed. 

Another interesting point with regard to figure 4(e) is the wide 
variation in phase angles at different flight speeds. Predicted trans- 
fer functions indicate that at low speeds (below a Mach number of 0.7) 
where ^ng a is negative, the phase angles approach -270° asymptotically 
with increasing frequency. Unpublished wind-tunnel data indicate that 

near a Mach number of 0.7 there is a transition in which Cn^ becomes 

-LL oa 

positive and consequently three of the coefficients in the numerator of 
the transfer function change sign. The result is an increase of l8o° 
in the high frequency phase lag. 


The frequency-response test points derived from flight data at a 
Mach number of 0-8l have been replotted in figure 5 which shows all six 
responses for the 35,000 foot altitude and the r/& r and p/6 a responses 
for 10,000 feet. These results are typical in indicating the degree of 
scatter usually encountered in the Fourier analysis of a particular 
flight record. Plotted as solid lines for comparison are predicted 
responses that have been calculated using estimates of the various 
stability derivatives presented in table II, which were obtained from 
reference k and also from wind-tunnel tests by the manufacturer . These 
calculations were made as shown in appendix C, using the exact linear 
fourth-order response equations. The agreement between measured and 
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predicted responses is generally good except, in some cases, at the 
extremes of the test frequency range. 

In addition to the pulse-type inputs, a step disturbance of the 
rudder was used in one instance at a Mach number of 0.8l in order to 
more clearly define the low frequency portions of the p/6 r and r/5 r 
frequency responses. Results of this analysis from 0.1 to 1.0 radians 
per second are plotted in figures 5(a) and 5(b) and are discussed in 
more detail in a later section. 

Also shown as dotted lines in figures 5(b), 5(c), 5(d), 5(g ) , and 
5(h) are responses computed from predicted transfer functions that have 
"been simplified as indicated in the following paragraphs. 


Theoretical Transfer Functions 


It is shown in appendix C that the characteristic equation A can 
he factored into the form 

A = D(D-X x ) (D-X 2 ) (D 2 +c x D+c 2 ) 

where X x and X 2 are the spiral and rolling roots, respectively, and 
where and c 2 are coefficients that define the oscillatory mode. 

By neglecting X 1 (which is usually very small) and by omitting other 
small terms that appear in the numerators of the Various response 
equations, three of the six responses may be reduced to the following 
simple forms: 


r _ B 3 D 

D ^4*0 j^D"t"C 2 


3_ = -B* 

5 r D 2 +CjP+c 2 

P_ = A 3 






( 1 ) 


It is also shown that by making additional assumptions as to rela- 
tive magnitudes, the coefficients X 2 , cj., and c 2 can be expressed as 
Lp, -(N r + Yp) , and Np*, respectively. Furthermore, since B 3 = N§ ’ 
and A 3 s B§ a ', equations (l) can be written as r 


r N 6r »D 

S r D^-(N r +Yp)D+Np' 


(2a) 
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3_ = - N S r 

5 r D 2 -(VY p ) D+ V 


(2b) 


!L W 

& a D-Lp 


(2c) 


It can be seen that in equations 2(a) and 2(b) the spiral and 
rolling modes are conpletely neglected, and yet, as shown in figures 
5(t>), 5(c), and 5(g ) , these simplified transfer functions yield responses 
that are almost identical to those obtained from the "exact" equations 
for frequencies greater than 1 radian per second. Similarly, the 
response computed from equation 2(c) closely matches the exact response 
(fig* 5(d)) over the frequency range shown except that it omits the 
small peak normally associated with the oscillatory mode. The spiral 
mode which has been neglected in all three simple equations appears to 
have no effect on the calculated airplane response except at frequencies 
well below 0.1 radian per second. 




Experimentally Determined Transfer Functions 


In the analysis of the flight data it was found that the frequency 
responses of r/8 r , |3/&r> and- p/§a could be successfully simulated by 
simple transfer functions of the same forms as equations (1). Solutions 
of these equations on the analogue computer, using final "best" values 
of the numerical coefficients with actual control motions as recorded 
i n flight, resulted in outputs that closely matched the measured time 
histories of r, 3, and p as shown in figure 3* This fact implies that 
the modes of motion that are neglected in each case have very little 
effect on the time response to a pulse-type input. 

By use of the coefficients that best describe the measured time 
histories, frequency responses were calculated for comparison with those 
derived directly from flight data. Examples of these calculations are 
shown by the dash-dot lines in figures 5(b), 5(d), 5(g), and 5(h). These 
curves, in general, match the experimental points closely for frequencies 
between 1 and 10 radians per second. 

The experimentally determined values of the coefficients X 2 , c x , 
and c 2 (using the notation of equations (l)) have been plotted in figure 6 
and are compared to predicted values of the same coefficients that were 
obtained by factoring predicted characteristic equations for several 
different Mach numbers. To give some idea of the errors involved in the 
assumptions of equations (2), predicted values of Lp, N r 4- Yp, and Np T 
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have also been plotted for comparison with the coefficients \ 2 , c lf and 
c 2 . It appears that, for this particular airplane, the assumptions are 
valid and that the simplified transfer functions form a logical basis 
for the evaluation of stability derivatives. The flight evaluated coef- 
ficients c x and c 2 have been transformed into the customary undamped 
natural frequency and damping ratio designations and are plotted in this 
form in figure 7* 


It is apparent from equations 2(a) and 2(b) that the same informa- 
tion can be obtained from either r/6 r or p/& r . Coefficients evaluated 
from each of these responses agreed favorably in most cases; however, 
because of indications that the yaw rate gyro possessed dynamic charac- 
teristics superior to those of the sideslip vane, only the yawing veloc- 
ity responses were used in the final calculations. 


The transfer functions of the three remaining responses p/8 r , r/& a , 
and p/b a were not amenable to simplification. However, it was found 
that the p/& r response could be matched satisfactorily by a transfer 
function of the type 


= 


5 r 


a 1 (D+a 2 ) (D+a 3 ) 
(D-A, 2 ) (d^dT^) 


which is the same form as developed in appendix C except that the spiral 
mode has been neglected. As written here, a x is identical to 
while a 2 and a 3 are complicated combinations of derivatives that can- 
not be readily simplified. Although this equation closely describes the 
measured time histories (fig. 3)> it was difficult to find unique values 
of the numerator coefficients. Changes in one of these could be compen- 
sated for by corresponding changes in the other two, and the values were 
not considered to be reliable enough for presentation. 

Definition of the r/b a and p/8 a responses required fourth-order 
transfer functions that include all three modes, and because of practical 
difficulties involved no attempt was made to evaluate the coefficients 
of these responses . 


Stability Derivatives 


In addition to the quantities Lp, N r + Yg, Ng T , and I^ 1 , that 
were determined as mentioned in the preceding paragraphs, the coefficient 
Lg 1 was evaluated from the time histories of rolling and yawing velocity 
as shown in reference 5* This method is briefly outlined in appendix D. 
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The quantities Ng, Lg, N^, and L§ a can then be calculated from 
Ng', Lg', N5 r ’ , and L§ a * by using the following expressions obtained 
from the relationships developed in appendix C: 


N 


3 


y-W 

1 _r x r z 


l q = 


= L 3 ,_r x N 3' 


1-r x r z 


Nc 


1 " r X r Z 



1_r x r z 


Because 
of Ls r ' were 
tion of Ng, r . 
ting L§ a . 


rv and are very small quantities, wind-tunnel estimates 
assumed to be sufficiently accurate to use in the calcula- 
The term r x ^5 ' was completely neglected in evalua- 

£L 


Finally, from the definitions given in the notation it is possible 

to evaluate the derivatives C 7 , C n , C 7 , C nc . , and C7 0 . 

*p n 3 ^3 lL b L b R 


The analysis methods used herein do not allow the separation of the 
damping term N r + Yp. As compared to N r , the term Yp is small and 
can generally be predicted accurately from wind-tunnel measurements . 
Therefore values of Cy^ given in table II were used in calculating 
Cn r from the quantity N r + Yp. 


The flight evaluated derivatives for both altitudes are plotted 
against Mach number in figures 8 and 9 * These are compared to the pre- 
dicted values listed in table II. Through the speed range of the test 
the predictions for both altitudes are essentially the same, except as 
noted in the plot of . 

The correlation between predicted derivatives and those evaluated 
from flight at 35,000 feet is generally good except, in some cases, at 
speeds near a Mach number of 1.0 where the predictions are apt to be 
inaccurate. Unpredicted variations with altitude are also apparent in 
the flight values of C^p, Cn r > and C n ^. At a Mach number of 0.8 the 
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value of Cn r for 10,000 feet is less than one half the value at 35,000 
feet. Flight data of reference 4 when expressed in this form show a 
similar trend. The plot of Cn 5 indicates the same tendency to a 

lesser degree, while the 10,000-foot value of Cjp (at M = 0.8) is some 
50 percent higher than the value at 35,000 feet. 

Values of Cjo determined in the present investigation agree favor- 
ably with wind-tunnel results, while those reported in reference 6 
(obtained from static flight tests of the same airplane) are much smaller 
in magnitude. It appears, however, that the results of reference 6 are 
subject to error because of the simplifying assumptions made. A more 
rigorous approach would have resulted in larger values of this derivative 


Examination of figure 9 shows the control effectiveness derivatives 

C n _ and C? ~ to have similar variations with increasing Mach number, 

k Og, 

and in each case the measured values are generally smaller than predicted 
Values of c obtained in the present investigation agree closely 


with those presented in reference 7 which again were evaluated from 
flight measurements of the same airplane. 


In this investigation there was no evidence of nonlinear variations 
of rolling or yawing- moment coefficients with p, r, or 3* This was 
concluded "because (l) the period and damping of the oscillations fol- 
lowing a control input were essentially constant in every case (no 
systematic variations with amplitude), and (2) the experimental time 
histories could he matched, in general, hy differential equations with 
constant coefficients. 


No conclusions are drawn as to nonlinear moment coefficient varia- 
tions with 5 a or 6 r because the magnitudes of the control inputs were 
not varied appreciably during the tests. They were small enough, however 
so that it could be assumed that the linear ranges were not exceeded. 


Aeroelastic Effects 


Although the present investigation was not conducted for the pur- 
pose of studying aeroelasticity, the test results do show the influence 
of structural deformation as indicated in the discussion that follows. 

Effects at high frequency .- The frequency-response measurements of 
p/6 a indicate a mode of motion at high frequencies that is not consis- 
tent with the rigid-airplane equations given in appendix C. According 
to these equations, the amplitude of this response should approach zero 
(at a slope of -1 on a logarithmic plot) as frequency increases, while 
the phase angle approaches -90° asymptotically. Examination of 
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figure 5(d), however, shows an increase in phase angle and a gradual 
decrease in amplitude attenuation at frequencies greater than 10 radians 
per second. These data indicate the presence of an additional high 
frequency mode which conceivably could correspond to the primary bending 
frequency of the airplane wing (about 8 cps) . A mode of this type can 
be included in the equations of motion, in general, by introducing 
additional degrees of freedom which will relate the mass and stiffness 
characteristics of the wing to actual motions of the airplane. From 
purely geometric considerations, bending of a swept-back wing is accom- 
panied by a change in angle of attack and, when the motions of the two 
wing panels are out of phase, there is a resultant rolling moment. Thus 
it is reasonable to expect wing bending or twist to have a noticeable 
effect on the rolling response of the airplane. 

Effects on flight evaluated stability derivatives .- Test values of 
^Zp shown in figure 8 indicate a variation not only with Mach number 
but also with altitude (dynamic pressure) which may be the result of 
inertia loading. 

When the airplane is accelerated in roll, inertia forces cause the 
two wing panels to bend in opposite directions. The resulting angle-of- 
attack variation produces a moment that tends to modify the rate of roll. 
This effect can be considered in the basic equations of motion by the 
introduction of dc^/dp, the variation of rolling-moment coefficient with 
rolling acceleration. Then letting 


qSb 8 Ct_ 

^ " ix dp 

equation (C 1 )-) may be written as 

j^(l-I^)D 2 -LpDjcp+(-r x D 2 -L r D)i|r-L33 = L§ a 6 a 

It should be noted that the inclusion of Lp is analogous to a change 
in the moment of inertia about the X axis. 

If each term is divided by (1 - Lp) then this equation is parallel 
in form to equation (Cl+), and the simplified transfer function for p/& 
becomes 

p _ k&a / l~lp 
& a " D-^Lp/l-Lp) 

Thus it appears that the two derivatives as evaluated from flight data 
are actually Cj p /(1 - Lp) and /(l " Lp) • For a swept-back wing 

L-p is a positive quantity that increases with dynamic pressure; so the 
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* measured values of the two derivatives should also increase with dynamic 
pressure. 


The test measurements of Cjp clearly show this trend, hut in the 
case of there was no appreciable variation with altitude. Theo- 

retical studies, however, and also wind-tunnel tests of flexible models 
(such as described in reference 8) have clearly indicated that, in the 
absence of inertia loads, the aileron effectiveness will in general 
decrease with increasing dynamic pressure, the primary contributing 
factor being the torsional flexibility of the wing. 


It has been shown in reference 8 that the derivative C^p will 
also vary with dynamic pressure. For straight wings there will generally 
be an increase in the negative value of C^p, while for highly swept- 
back wings there will be a decrease. Presumably there is an intermediate 
sweep angle at which there is essentially no variation in Cj with 
dynamic pressure. 


In order to substantiate this hypothesis regarding the effects of 
elasticity, numerical calculations were made using the method of refer- 
ence 9 to estimate the variations with dynamic pressure of the quanti- 
ties C^p, , and Lp. These calculations were based on information 

supplied by the manufacturer regarding the mass distribution and stiff- 
ness of the airplane wing. While not agreeing quantitatively with the 
measured variations, the results did show the same trends with dynamic 
pressure that were observed in the flight evaluated derivatives. More 
specifically, the calculations indicated that Cj decreases much more 

°a 


rapidly with increasing dynamic pressure than does Cj . Furthermore, 
the variation of (l - Lp) was approximately the same as that of • 

This indicates that an increase in the effective derivative C^pAl - Lp) 
with dynamic pressure can be expected even though there is little varia- 
tion in /(I - L]5) • 


Similar arguments may be advanced to explain the altitude variations 
that were noted in C nr and C n ^ . These variations could be the result 

of distortions of the fuselage and tail due to inertia and aerodynamic 
loadings . 


Response at Low Frequency 


As mentioned in an earlier section, responses to a rudder step 
disturbance were recorded in flight in order to check mainly the low 
frequency portion of the r/5 r frequency response (0 to 1.0 radians per 
second). As shown in figure 5(b), the results of this analysis verify 
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the prediction of a sharp attenuation in amplitude at a frequency of 0.4 
radians per second. A step input was used in this particular case because 
it gives much more accurate results at low frequencies than a pulse input. 
The question of control inputs is discussed more fully in appendix B. 

Although not shown in figure 5(b) , the zero frequency (steady-state) 
amplitude of r/b r was measured to be 1.44 as contrasted to a value of 

1 69 .0 predicted from the exact equation. Similarly, the steady-state 
amplitude of p/b r was measured as 1.75j while the theoretical equations 
indicate that this value should be zero. These discrepancies can be 
partially explained in view of the following discussion. 

The linear equations of motion used herein are valid only for small 
angular displacements. When the step disturbance was applied in flight, 
an angle of bank of approximately 70° was reached before the motions of 
the airplane became steady. In the side force equation (equation (C6)) 
the term K^cp is actually a linear approximation of K 2 sin cp and 
obviously is valid only for small bank angles. To determine the effect 
of this nonlinearity on the predicted airplane responses, time histories 
of p, r, and 9 were calculated on an analogue computer for a rudder 
step input of the same magnitude as applied in flight. Predicted values 
of the various stability derivatives were used and solutions were obtained 
first with the linear approximation K x cp and next with the nonlinear 
K-l sin Cp. Results obtained from the nonlinear equation indicated that 
the rolling velocity response would reach a finite steady-state magnitude 
as seen in the flight record. In the case of the yawing velocity 
response, however, the effect of the nonlinearity was such as to reduce 
the steady-state value somewhat but not nearly enough to account for the 
extreme difference between flight measurements and predictions. 

The relationship of spiral damping to the zero frequency amplitude 
offers another possible explanation for this discrepancy. Although the 
basic analysis methods used in this report give no information concerning 
the spiral mode, predictions indicated that at a Mach number of 0.8 at 

35.000 feet the spiral root was -0.00070 which leads to the previously 
mentioned steady-state magnitude of 169 .O for r/& r . The subsidence of 
this mode was clearly measurable, however, in the yawing velocity response 
to the rudder step input. A line was drawn through the center of the 
free oscillations to represent the spiral mode and the root was measured 
as -0.07 instead of -0.00070 as predicted. 

This measurement was verified by examining recorded time histories 
of yawing velocity responses to aileron pulse inputs (e.g., fig. 3(b)). 

On eight flight records taken at Mach numbers between 0.5 and 0.9 the 
spiral root was found to vary from -0.06 to -0.09. It can be seen that 
the discrepancy in spiral subsidence is of the same order of magnitude 
as the discrepancy in steady-state magnitudes of r/B r . 


CONFIDENTIAL 


NACA RM A52I17 


CONFIDENTIAL 


19 


In the characteristic equation the coefficient C Q is generally 
very small as compared, to C-,^ and thus the theoretical spiral root can 
he expressed quite accurately as simply 

(3) 


which is approximately the same as 




K x (L r K 3 -N r L 3 ) 


(M 


All the quantities on the right side of this expression except L r have 
been determined experimentally and found to agree reasonably well with 
predicted values. Even though the term L r Np - N r Lp represents a small 
difference of two large numbers, it is inconceivable that the errors in 
Np, N r , or Lp can be of sufficient magnitude to account for the large 
deviation in the spiral root. Therefore it appears that this discrepancy 
must be at least partially due to an erroneous estimate of L r . By 
rearranging terms in equation (4) 


L r 


^l-kp N r Lp 
■ Ki + H p 


Using experimental values of all quantities on the right-hand side for 
the case of 0.8 Mach number at 35>000 feet results in a value of -4.2 
for L r as contrasted to the original prediction of 0.844. Further 
study showed that the resulting value of L r was not appreciably changed 
by using more exact expressions in place of equations (3) and (4). Then 
assuming that all experimentally determined values on the right-hand side 
of equation (5) are correct, Ci r must have a value of -0.537 rather than 
0.108 as predicted. 


Effects of Minor Derivatives, Product of 
Inertia, and Flight -Path Angle 


The results of this report indicate that a knowledge of C n p> 

C nr , and Cy-p is sufficient to define the oscillatory mode, that the 
rolling mode is primarily a function of Cj , and that the spiral mode 
is defined by C^, C n(3 , C nr , C^, and C^. Therefore it would appear 

that C n has little effect on the dynamic behavior of the test air- 
plane and that Ci T influences only the spiral mode. Because deriva- 
tives such as C n p and Cj r are usually difficult to estimate, the 
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question may arise as to whether large changes in either of these deriva- 
tives will appreciably alter the airplane stability. If erroneous esti- 
mates of and Ci r had been used in this report, would the simplifi- 

cations made in appendix C still be valid? Similar questions may arise 
regarding the effects of product of inertia and flight -path angle. For 
this reason table III has been prepared. This table shows the manner in 
which radical changes in each of the above-mentioned quantities will 
effect the roots of the characteristic equation. In this example, theo- 
retical data for a Mach number of 0.5 and an altitude of 35*000 feet were 
used. Sample calculations have shown that at this speed and altitude the 
effects to be considered are more extreme than for any of the other flight 
conditions covered in this investigation. 

Examination of table III shows trends that have been verified in 
many other investigations of this type. Extremely large changes in 
Cnp result in variations of all four of the roots although only the 
oscillatory damping is greatly affected. On the other hand changes in 
Ci , as previously intimated, cause large variations in only the spiral 
damping. Variations in 1^ appear to influence all the roots except 
the spiral root, with the largest effect on the oscillatory damping. 
Finally the introduction of a flight-path angle into the equations changes 
the spiral root radically, but even for an angle of -90° there is very 
little effect on the other roots. 

Looking at table III from a different point of view, it can be seen 
that the oscillatory damping is influenced to some extent by each quan- 
tity considered and thus is often difficult to predict accurately. This 
obviously means that the expression N r + Yn may not always be adequate 
in defining the oscillatory damping. 


CONCLUDING REMARKS 


A flight investigation has been performed on a 35 ° swept -wing fighter - 
type airplane in which dynamic lateral- and directional -response charac- 
teristics were measured. Transient -type responses to rudder and aileron 
disturbances were recorded at altitudes of 10,000 and 35,000 feet in the 
Mach number ranges of 0.50 to 0.8l and 0.50 to 1.04, respectively. From 
the results of the analysis of these data, the following statements can 
be made. 

Airplane responses in yawing velocity and side-slip angle due to 
rudder disturbances can be represented by second-order transfer functions 
that are related solely to the oscillatory mode. Simple first-order 
equations adequately define the rolling velocity response to an aileron 
input. It was found that these equations would closely define an entire 
measured time history and also describe the corresponding frequency 
response through the range of 1 to 10 radians per second. 
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Fourth-order transfer functions calculated from the basic equations 
of motion using wind-tunnel and theoretical estimates of the various 
stability derivatives can be simplified by neglecting small quantities 
and by making approximate cancellations until they are of the same form 
as those evaluated from flight data. Furthermore, it was possible to 
express the coefficients of these transfer functions in terms of indi- 
vidual stability derivatives. Frequency responses computed from these 
simplified equations were almost identical (between 1 and 10 radians per 
second) to those conputed from the exact fourth-order transfer functions, 
and when compared with experimental results there was generally good 
agreement. Thus it is concluded that the simplified transfer functions 
form a reliable basis not only for estimating airplane responses but 
also for the flight evaluation of stability derivatives. The methods 
used here are felt to be sufficiently general to apply to any conven- 
tional airplane, with some reservation regarding the accuracy of the 
evaluation of Cn r • 

Experimental values of the derivatives C^p, Ci p, C^p, Cn r > 
and Ci~ compared favorably with predictions, based on theory and wind- 
tunnel measurements , at Mach numbers below 0-95, while at higher speeds, 
where predictions are questionable, there was some deviation. There were 
also notable discrepancies in flight values of Cjp and C nr , obtained at 
the 10,000 foot altitude, which were attributed to structural deforma- 
tions resulting from aerodynamic and inertia loads. 

When the. Fourier analysis for p/& a was expended ^ 0 frequencies 
beyond 10 radians per second, the frequency response showed evidence of 
aeroelastic deformation which appeared as an additional Aode of motion 
not consistent with rigid airplane theory. 

As the frequency approaches zero the spiral mode becomes the pre- 
dominant factor in the airplane response. The spiral root was measured 
as -0.07 (at a Mach number of 0.8) which is many times greater than 
predicted. To satisfactorily account for this large discrepancy would 
require a negative value for C Zr which is contrary to theoretical 
estimates . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, California 
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APPENDIX A 
TRANSFER OF AXES 


The equations of motion normally used in airplane dynamics are 
"based on a system of axes fixed in the airplane in which the X axis is 
the intersection of the plane of symmetry and a plane perpendicular to 
the plane of symmetry that contains the relative wind vector. These are 
normally referred to as stability or flight -path axes. The angular dis- 
placement between the X axis and the reference axis of the airplane is 
equal to the angle of attack. Since recording instruments are generally 
alined with the reference axis, measurements of angular displacements 
and rates must be corrected to conform to stability axes notation as 
indicated in the following sketch taken from reference 5- Here <p and 
are vector components of the resultant rotation of the airplane and 
the subscript 1 refers to the reference or body axes. 



From the sketch it can be seen that 

(p = cp x cos a + \|r x sin a 
^ ^ cos a - cp x sin a 

Sideslip angles can be transformed by the relation 

3 = cos a 

For most purposes these conversions need be made only when the 
angle of attack is large. In this investigation it was found that the 
corrections to rolling velocity responses could be neglected in «i 1 cases 
because, for either type of input, the response in yaw is small compared 
to that in roll. In the case of the yawing velocity records, however. 
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corrections were necessary. When a rudder input was used this correction 
was negligible at those speeds where angles of attack were less than 2°, 
but in the case of the yawing responses to aileron inputs the conversion 
had to be made at all speeds . The correction to 3 was neglected in 
every case. 
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APPENDIX B 
CONTROL INPUTS 


When frequency responses are to be calculated from transient records, 
care should be given to the choice of a suitable forcing function. The 
frequency range through which accurate transformations can be obtained 
is definitely limited by the shape of the control input. Theoretically, 
a pure impulse (zero time duration) is the most desirable input for all 
purposes because it gives uniform excitation to the entire frequency 
spectrum. The transform of a step input, on the other hand, has a mag- 
nitude that varies inversely with frequency and thus gives infinite 
excitation to the zero frequency component at the expense of the higher 
frequencies . 

The nearest physical approach to a pure impulse is an input that is 
roughly triangular in shape as shown in the sketch: 



Letting a equal the slope and T equal the time base of the triangle, 
the Fourier transformation of this input can be obtained from the relation 


S(ioj) 


= / S(t)e“^ u ’^dt 


^ o 

Integration results in a transformation with the following real and 
imaginary parts: 



(Bl) 

(B2) 
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The magnitude of the transformation is then 



It can be seen that | S | is periodic and is zero when w = kn/l, 

8k /T, At these frequencies the transform of the response to this 

input would also be zero, and thus the ratio of output to input would be 
indeterminate. A reduction in T would increase the period and reduce 
the number of indeterminate points . This is shown in figure 10 where 
the transform magnitudes of two triangular pulses are plotted. One has 
T = 1 second and a = while the other has T = l/2 and a = 16. The 
areas under the two triangles are equal so that their transforms have 
equal magnitudes at zero frequency. Reducing T from 1 second to 1/2 
second doubles the period and moves the first indeterminate point from 
a frequency of 4 k radians per second to 8k. For purposes of compari- 
son, transformations of a unit step and a unit impulse are also shown. 

As T is further reduced, the magnitude of the triangular pulse 
more closely approaches the constant value that is characteristic of the 
pure impulse. To gain full advantage from the smaller T, the slope 
must be increased to maintain the same area under the pulse. A practical 
limitation is fixed by the maximum rate at which a control surface can 
be moved, and any further reduction in T results in smaller over-all 
magnitudes. The most desirable input, therefore, is a compromise between 
large area and short time duration. 

From figure 10 it would appear that a pulse-type input is well 
suited for determining low frequency characteristics. However, the 
following explanation will show that this is not true. 

Generally it is impossible to return a control surface precisely to 
its initial position after application of a pulse input. Even if a chain 
stop or other device is used there is still apt to be a small residual 
deflection after time T as shown in the following sketch: 
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If 5 2 is exactly zero, then as w approaches zero (from equations 
(Bl) and (B2) ) 


Now if 5 2 is finite, the real and imaginary parts of the transform of 
the entire input are 


thus the zero frequency magnitude is infinitely large regardless of how 
small 6 2 may he. Therefore, even though & 2 appears to be zero on a 
flight record (i.e., 8 2 is less than the least count of the recording 

instrument) there is still the possibility of an infinite error at zero 
frequency. A step input is not subject to these large low frequency 
errors; an error of 1 percent in the reading of the step deflection 
merely means an error of 1 percent in the transformation. 





In this case as w — > 0, 


I 


00 
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APPENDIX C 

PREDICTED AIRPLANE RESPONSES 
EQUATIONS OF MOTION 


The three equations that define the lateral and directional motions 
of an airplane with respect to stability axes such as developed in refer- 
ence 10 may be written as: 

(l x D 2 - qSb C Zp ~D)q> +(dxzD 2 - qSb C Zr ~ d) * - qSb C^p = qSb C Z& S 

(Cl) 

(-I XZ D 2 * <lSb C^ |^d) cp + ^D 2 - qSb C nr |^d) * - qSb C np p = qSb C n& & 

(C2) 

(-W cos 7)9 + (mVD - W sin 7)* + (mVD - qSCy p )p = qSCygb (C3) 


By dividing equation (Cl) by I x , equation (C2) by I z , equation (C3) by 
mV, and by introducing new symbols, the three equations can be written 
in the more convenient form that follows: 


(D 2 - LpD)cp + (-r x D 2 
(-r z D 2 - NpD)cp + (D 2 
-K x 9 + (D - K 2 )^ 


- L r D)V - Lpp = L&5 

(C4) 

- N r D)V - NpP = N 5 5 

(C5) 

+ (D - Y p )p = Y & & 

(c6) 


CHARACTERISTIC EQUATION 


The characteristic equation A is formed by expanding the major 
determinant to give 

A = D(C 4 D 4 + C 3 D 3 + CgD 2 + C X D + C 0 ) (C7) 
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where 


c 4 = 1 - r X r Z 


C 3 = -Lp - N r - Yp(l - r x r z ) - r x N p - rfr T 

C 2 = (N r + r^LrJYp + (Lp + r x Np)Yp + (LpN r - NpL r ) + (Np + r z Lp) 

C x = -(LpN r - NpL r )Yp + (LpNp - NpLp) - K x (Lp + r x Np) - K 2 (Np + r z Lp) 

C Q = -K x (L r Np - NpLp) - K 2 (LpNp - NpLp) 

These coefficients can be further simplified, by making the following 
substitutions : 

Let 


Then the .coefficients of the characteristic equation are finally expressed 
as 


C 3 = -Lp* - N r » - Yp(l - r x r z ) 

C 2 = (N r * + Lp*)Yp + (LpN r - NpL r ) + Np* 

. Cl = -(LpN r - NpL r )Yp + (LpNp - NpLp) - K x Lp» - K^p* 
C 0 = -Ki(L r Np - N r Lp) - K 2 (LpNp - NpLp) 

In factored form, equation (C7) is 


where \ x and \ 2 are designated as the spiral and rolling roots, respec- 
tively, and where X 3 and X, 4 are a complex pair (a ± iu> x ) that describe 
the oscillatory mode. 

For convenience in this investigation A has been expressed as 



N r * — Np + r z Lp 
Np’ = Np + r z Lp 


C 4 = 1 - r x r z 


A = D(D - X x ) (D - X 2 ) (D - Xg) (D - \ ) 


A = D(D - Aa) (D - X 2 ) (D 2 + c x D + c 2 ) 
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where 


c x = ~(\ 3 + \ ) = -20 
C 2 = ^3^4 = a 2 + Wl 2 

Here c-l and c 2 are real coefficients that define the damping and period 
of the oscillatory motion. 


The quadratic term may also he written in the form 


* fl + 2{i + 

11 V “n "n/ 


where 


and 


»n " ^ 

C 1 


SIMPLIFICATIONS OF THE CHARACTERISTIC EQUATION 


In the special case when the flight-path angle 7 is zero and when 
the product of inertia is very small, it is often possible to neglect 
Co and, by neglecting other small terms, the characteristic equation 
may be written as 

A = D 2 Jd 3 - (Lp + N r + Yp)D 2 + (LpYp + LpN r + Np)D - LpNp 
The cubic term can be factored exactly so that 


A = D 2 (D - Lp) 




This form of the characteristic equation considers only the oscillatory 
and rolling modes. It enables the coefficients C 2 , C3 to t>e 

expressed directly in terms of aerodynamic derivatives or simple com- 
binations thereof . 


Even when the product of inertia is significant the characteristic 
equation may be factored approximately into the comparable simple form 


A = D 2 (D - Lp) 



(N r + Yp)D + Np* 


While the factorization is not exact, it is nevertheless justifiable 


in many cases . 
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TRANSFER FUNCTIONS 


From the three equations of motion the airplane responses in cp , 
if, and 3 can he readily calculated. In the following equations, 6 

refers to either an aileron or a rudder disturbance . 

% 

cp A 3 D 3 + AgD 2 + A-lD + A 0 
6 = A 


where 


A 3 = l 5 + r X N S 

A 2 = -Lg(Yp + N r ) + Ng(L r - r-^Yp) + YgL p * 

Ai = Lg(N r Y p + N p ) - Ng(L p + L r Y p ) + Yg(L r N p - N r L p ) 
A 0 = K 2 (NgL p - N p Lg) 


* B3D 3 + B2D 2 + BtD + B n 
S A 

where 

B 3 = N S + r Z L 6 

B 2 = -Ng(Lp+Y p ) + Lg(N p - r z Y p ) + YgN p » 

B z = %S Y 3 “ L 5 N p Y 3 + Y 5( L P N p - ^Lp) 

B 0 = K 1 (L 5 N p - NgL p ) 


P D(E3D 3 + E^p 2 + EjD + E 0 ) 

6 " A 

where 

E 3 = y 6 ^ B " r X r Z^ 

E 2 = "^(Np 1 + Ep'^ " B S r Z - E S 

= -^(Ep^r “ ^pl'r) " Bg(Np - r z K 2 - K 1 ) + ^(r-^K^ + Lp + K 2 ) 
E 0 = KxCNgLp - LgN r ) + K 2 (LgN p - NgLp) 
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Simplifications can also be made in these expressions by neglecting 
small quantities ; however, this can be shown more clearly in the numeri 
cal example that follows. 


NUMERICAL EXAMPLE AT M = 0.8 


Using values of stability derivatives shown in table II, and with 
7 = 0, predicted responses for p = Dcp, r = D^, and 3 f° r both aileron 
and rudder inputs have been calculated and found to be 

A = D(D 4 + 3-652D 3 + 15.16D 2 + 41.26D + 0.0289) 

= D(D + 0.00070) (D + 3*078) (D 2 + 0.573D + lS-^O) 

_r_ _ -7.60D(D + 3.091) (D 2 + 0.0270D + 0.208) 

5 r A 

_p_ 5.16D 2 (D + 4.U36) (D - 5.210) 

6 r A 


3_ = 0.0339D(D + 3-053) (D + 225.3) (p - 0.00703) 
S r A 

r_ = 0.699D(D + 3.978) (D s - 1.798D + 7.398) 

®a A 


_p_ = 36.4 d 2 (d 2 + 0.65SD + 13.68) 

5 a ~ A 


_3_ = 0.0008 d(D + 0.990 ) (D - 1.09*0 (d - 870) 

s a " A 

By neglecting small terms, r/b r can be expressed as 


r = -7.60D 3 (D + 3-091) 

5 r D 2 (D + 3.078) (D2 + 0.573D + 13.40) 


1 
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and then "by an approximate cancellation this reduces to 


r _ -7.60D 

6 r D 2 + 0.573D + 13.40 

Similarly, 3/& r can be simplified by neglecting small terms so that 


J3_ = 7.64d 2 (d + 3-053) 7.64 

S r D2(D + 3.078) (D2 + 0.57 3D + 13.40) ~ D 2 + 0.573D + 13.40 

It can be seen that this expression for p/5 r is practically indentical 
(with opposite sign) to the integral of the simplified equation for 
r/& r . It is also possible to simplify p/& a as follows: 

p _ 36.4 e^(d 2 + 0.655D + 13.68) ^ 36.4 

6 a D 2 (D + 3.078) (D2 + 0.573D + 13*40) ~ D + 3-078 

Similar simplifications have been made for other Mach numbers and 
found to be equally valid. 


/' 

k 
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APPENDIX D 

FLIGHT EVALUATION OF C 7 


From time histories of the free oscillatory responses of p and r 
it is possible to evaluate the derivative C 2 provided the derivatives 

C n p and are known. The method employed is explained in some detail 

in reference 5 and is briefly summarized here. 


When the three equations of motion are set equal to zero and written 
in determinant form, expansion of the appropriate minor determinants 
about the third row yields cofactors of 9 and V . 


The cofactor of 9, Qp is 


Qp = (r x D2 + L r D)N p +.(D2 - NjDjLp 


= D 


Lp'D + 


(L r N p - N r Lp) 


The quantity LpNp - NpLp is generally very small and can be closely 

approximated as 


Cep = Lp’D 2 


The cofactor of Cty is 

Cty = -(D 2 - LpD)Np - (r z D 2 + N p D)L p 
= -dL’D - (LpNp - NpLp) 

v L 

In this case the term NpLp is small compared to IpNg and the expres- 
sion for Cty can be approximated as 

C\|c = -DNp*(D - Lp) 

The ratio of the two cofactors is then 

C9 _ Lp r D 
C t -Np«(D - Lp) 

When the complex root = a + iu 1 is substituted for the operator D, 
this expression is the ratio of the free oscillatory responses of 9 and 
t at any time t. The ratio of p to r is obviously the same and can 
be expressed as 
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P 


-Lg* (c + 


r 


N p * (o’ + iw x - L p ) 


The actual magnitude of this ratio is 


P 


Lg* t/cf 2 + 


r 


N p * */( a “ Lp) 2 +U)! 2 


In this form |p/r| is the ratio at any time t of 'the amplitudes of 
the envelopes that enclose the oscillatory motions of p and r; a is the 
rate of damping of the envelope; and is the natural frequency of 

oscillation. When a is very small as compared to it can be omit- 

ted; thus 


If Np* and Ip are known, it is then possible to evaluate Lp* from 
measured time histories of p and r. 


P 


Lp’gj], 


Lp 


r 



Np'/l + (Lp/ W i ) 2 
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TABLE I.- PHYSICAL CHARACTERISTICS OF TEST AIRPLANE 


Wing 


Total area 

Span 

Aspect ratio 

Taper ratio 

Mean aerodynamic chord 

Dihedral ► . * . . 

Sweepback of quarter- chord line 

Aerodynamic and geometric twist 

Root 'airfoil section (pormal to quarter- 

chord line) 

Tip airfoil section (normal to quarter- 
chord line) 


287.9 sq ft 
. . 37.1 ft 
. • • 4.79 
. . . 0.51 
. 97.03 in. 
.... 3° 
. . 35°l4* 
. . . . 2 ° 


NACA 0012-64 (modified) 
NACA 0011-64 (modified) 


Ailerons 


Area, each 18.6 sq ft 

Span 9*18 ft 

Chord, average 2.03 ft 

Deflection, maximum 14° up, l4° down 

Inboard end at 51*6 $b/2 

Vertical tail 

Area, total 34.4 sq ft 

Span 7-5 ft 

Aspect ratio 1.74 

Taper ratio O.36 

Sweepback of quarter-chord line 35°00’ 

Rudder 

Area ' 8.1 sq ft 

Span 6.6 ft 

Chord, average 1.23 ft 

Deflection, maximum 24.8° right, 25° left 

Average weight for calculations 12,800 lb 

Moment of inertia about X axis . 7>245 slug ft 2 

Moment of inertia about Z axis 23,190 slug ft 2 

Inclination of principal longitudinal axis 

with respect to fuselage reference axis -2.5° 
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TABLE II.- PARAMETERS USED IN ESTIMATING AIRPLANE RESPONSES 



Altitude, 35,000 feet 

Alt: 

itude y 

10,000 

feet 

M 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

0.5 

0.6 

0.7 

0.8 

V 

486 

583 

681 

778 

875 

972 

538 

646 

754 

862 

<1 

87.O 

125.2 

170.6 

222.5 

281.8 

347.5 

25^.7 

366.7 

499.0 

652.0 

cl 

.51 

.36 

.26 

.16 

• 13 

.12 

.17 

.12 

.09 

.07 

a 

7.2 

5.3 

3-8 

2.8 

2.4 

1.2 

2.1 

1.3 

.8 

.4 

Ixz 

-1297 

-773 

-359 

-83 

28 

386 

124 

331 

469 

580 

Cl B 

-.1025 

-.0857 

-.0773 

-.0741 

-.0721 

-.0768 

-.0627 

-.0573 

-.0563 

-.0571 

c ng 

.1100 

.1146 

.1199 

.1273 

.1366 

.1467 

.1105 

.1144 

.1200 

.1270 

c 4 

-.690 

-.701 

-.715 

-.733 

-.757 

-.782 

-.690 

-.701 

-.716 

-.735 

c 2p 

-.360 

-.367 

-.375 

-.385 

-.399 

-.414 

-.358 

—.366 

-375 

-.386 

C -P 

-.0328 

-.0225 

-.0160 

-.0120 

-.0092 

-.0068 

-.0113 

-.0076 

-.0055 

-.0044 

Cll 

• 157 

.130 

.116 

.108 

.106 

.104 

.091 

.082 

.078 

.076 

Cn r 

-.1820 

-.1852 

-.1896 

-.1970 

-.2065 

-.2170 

-.1817 

-.1866 

-.1922 

-.1984 

C Z 8 r 

.0077 

.0102 

.0138 

• 0155 

.0160 

.0183 

.0130 

.0160 

.0180 

.0200 

Cn§ 

-.0730 

-.07?8 

-.0725 

-.0742 

-.0736 

-.0582 

-.0729 

-.0728 

-.0733 

-.07*t2 

CY 5r. 

.160 

.160 

.160 

• l60 

.160 

.160 

.160 

.160 

.160 

.160 

r ° r 

.112 

.114 

.113 

.111 

.088 

.043 

.110 

.110 

.110 

.110 

u a 

Cn 6 a 

-.0050 

.0020 

.0059 

.0081 

.0095 

.0105 

.0050 

.0095 

.0110 

.0120 

C Y6 a 

.004 

.004 

.004 

.004 

.004 

.004 

.004 

.004 

.004 

.004 

Xa 

-.00192 

-.00113 

-.00076 

-.00070 

-.00027 

-.00077 

-.00144 

-.00153 

-.00146 

-.00123 

A >2 

-1.809 

-2.203 

- 2.667 

-3.078 

-3.581 

-4.168 

-4.77 

-5.87 

-7.01 

-8.17 

C 1 

. 378 

.438 

.497 

• 573 

.679 

.7to 

.910 

1.119 

1.322 

1.524 

C 2 

5.24 

7.25 

9.91 

13. *«> 

18.03 

.23.45 

12.75 

19-48 

28.35 

38.35 
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TABLE III.- EFFECTS OF Np, L r , lyz, AND 7 ON 
COEFFICIENTS OF THE CHARACTERISTIC EQUATION 

[A = (D - X 2 ) (D - X 2 ) (D a + c x D + c 2 )] 


(a) Effect of variations in Cnp 

% 

Spiral 

root 

Rolling 

root 

x 2 

Oscillatory 

mode 

c i 

c 2 

- 0.10 

-.05 

0 

.05 

-0.00168 

-.00192 

-.00206 

-.00223 

- 1.951 

- 1.809 

- 1.731 

- 1.646 

0.235 

.378 

.455 

.541 

5.55 

5.24 

5-10 

4.96 

(b) Effect of variations in Cj r 

L r 

Xi 

X 2 

c i 

C 2 

1.688 

.844 

0 

-.844 

0.0216 

-.0019 

-.0262 

-.0520 

- 1.818 

-I.809 

- 1 . 8 o 4 

- 1-795 

0.393 

.378 

.359 

.342 

5-30 

5.24 

5.18 

5.12 

(c) Effect of variations in Ixz 

I XZ 


X 2 

Cl 

c 2 

-2594 

-1297 

0 

-0.00191 

-.00192 

-.00193 

-1.680 

- 1.809 

- 1.956 

0.608 

.378 

.176 

5 - 86 * 

5.24 

4.78 

(d) Effect of flight-path angle 

7 



C 1 

C 2 

1 1 

VO uo 

000 
0 0 

-0.00192 

-.0313 

-.0595 

-I.809 

-1.812 

-I.816 

0.378 

.356 

.314 

5.24 

5-24 

5.23 
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Figure 1.— Photograph of test airplane. 
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Rudder deflection, Yawing velocity, Rolling velocity, p, radians/sec Sideslip angle, 

S r , deg r, radians/sec ft, deg 
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kl 



.2 


0 


-.2 




Time, t, sec 

(a) Rudder pulse -type input. 


Figure 3. - Sample flight records of yawing velocity, rolling velocity, 
and sideslip angle at a Mach number of 0.8/ . 
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(b) Aileron pulse - type input. 
Figure 3- Concluded. 
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Frequency, to, radians/sec 

(a) Rolling velocity response to rudder input. 

Figure 4.- Lateral - directional frequency responses at various flight 
Mach numbers at 35,000 feet. 
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Frequency, cu, radians/sec 
(b) Yawing velocity response to rudder input . 
Figure 4 - Continued . 
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Frequency, < 0 , radians/sec 
(c) Sideslip angle response to rudder input . 
Figure 4 -Continued. 
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Frequency, co, radians/sec 
(d) Rolling velocity response to aileron input. 

Figure 4 - Continued. 
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Figure 4. - Continued. 
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Frequency, oj, radians/sec 
(f) Sideslip angle response to aileron input. 
Figure 4.- Concluded . 
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(a) Rolling velocity response to rudder input , 
altitude, 35,000 feet. 

Figure 5 - Comparison of experimental and predicted lateral- 
directiona! frequency responses at a Mach number of 0.8/ . 
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Frequency, <u, radtans/sec 

(b) Yawing velocity response to rudder input, altitude, 35,000 feet 

Figure 5 - Continued. 
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Frequency, cu, radians/sec 
(c) Sideslip angle response to rudder input, altitude, 35,000 feet. 

Figure 5- Continued. 
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O Experimental data - pulse input 

Predicted response from exact equation 

Predicted response from simplified equation 

Response calculated from experimentally 

determined transfer functions 



(d) Rolling velocity response to aileron input; 
altitude, 35,000 feet. 


Figure 5- Continued . 
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O.l 1.0 10 20 


Frequency, < 0 , radians/sec 



(e) Yawing velocity response to aileron input ; 
altitude, 35,000 feet. 


Figure 5 - Continued. 
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Frequency, a>, radians/sec 


(f) Sideslip angle response to aileron input ; altitude, 35,000 feet. 
Figure 5- Continued 
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o Experimental data- pulse input 

Predicted response from exact equation 

Predicted response from simplified 

equation 


Response calculated from experimentally 

determined transfer functions 



Frequency, <u, radians /sec 

(g) Y a wing velocity response to rudder input; altitude, 10,000 feet. 
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Figure 5- Continued . 


VJI 

VJl 


KACA EM A52I17 CONFIDENTIAL 


5 6 


CONFIDENTIAL 


NACA EM A 52117 


<b 



0 Experimental data -pulse input 
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equation 

Predicted response from simplified 

equation 

Response calculated from experimentally 
determined transfer functions 
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(h) Rolling velocity response to aileron input ; 

10,000 feet. 
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Figure 5.- Concluded . 
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Figure 6- Variation with Mach number of the coefficients of 
the characteristic equation A £ , c /} and C e . 
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Figure 6. - Concluded . 
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Figure 7- Variation of undamped natural frequency and damping 

ratio with Mach number. 
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Figure 9- Variation of the derivatives C. and C n with 
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Figure to- Fourier transformation magnitudes of a step, an impulse, and two triangular shaped inputs . 
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